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AN ANALYSIS OF THE TRANSITION OF A HELICOPTER FROM 
HOVERING TO STEADY AUTOROTATIVE VERTICAL DESCENT 
By A. A. Nikolsky and Edward Seckel 

SUMMARY 


An analytical study is presented of the transition from hovering 
flight (power on) to steady vertical descent (in autorotation) following 
power failure while hovering. The effects of hinging the "blades, of "blade 
moment of inertia, and of rapidity of pitch reduction after power failure 
are investigated. - 

The results indicate that the effect of "blade flapping is negligible 
as far as the establishment of steady autorotation is concerned. From 
the standpoint of avoiding excessive "blade stalling during the transition, 
it is desirable that the "blade pitch "be reduced as rapidly as possible 
after power failure, and that the blade moment of inertia be large. 

Sample calculations and graphs of computed transitions for a typical 
helicopter are presented. 


INTRODUCTION 


This report is the second phase of a broad program of study of the 
transient motions of helicopters in autorotative flight . The first 
phase (reference l) dealt with the steady-state condition of autorotative 
vertical descent. This report is concerned with the transition from the 
steady condition of hovering to the steady autorotative descent. 

It has long been realized that the ability of a helicopter to 
establish steady autorotation following power failure would be an 
important safety feature. Although, some helicopters have successfully 
demonstrated such an ability, there has been no analytical method avail- 
able for assessing the influence of the design variables on the eaBe and 
safety of performing such a maneuver. 

Such a method, based on numerical integration of the equations of 
motion, 1 b presented in this report. 
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Validity of the analysis presented here is limited to cases where 
Hade stalling can he neglected throughout the maneuver, and therefore 
where the rotor is in no danger of stopping <. Nevertheless the investi- 
gation is sufficient to indicate the effects of certain important 
variables . 

This work was conducted at Princeton University under the sponsor- 
ship and with the financial assistance of the National Advisory Committee 
for Aeronautics » 
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■SYMBOLS 

Physical Quantities 


gross weight of helicopter, pounds 
number of blades per rotor 


blade radius, feet 

radial distance to blade element, feet 
blade -Section chord, feet 


equivalent blade chord, feet 



rotor solidity ratio 




blade-section pitch angle from zero lift, 
radians unless otherwise stated 

mass moment of inertia of blade about flapping 
hinge, slug-feet^ 

mass density of air, slugs per cubic foot 
time, seconds 

acceleration of gravity ( 32*2 ft/sec^) 


g 
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Y 

V v 

fl 


X 


Velocities 

true airspeed of helicopter along flight path, 
feet per second 

vertical component of Y (positive down, as 
in descent 

rotor angular velocity, radians per second 

induced velocity at rotor (always positive), 
feet per second 

inflow ratio (assuming v constant over the 



c d 


o 


Blade-Element Aerodynamic Characteristics 

section profile-drag coefficient 


e,- 


1 


&0,S 0 ' > & 1> & 2 


Ac d o ' A5 0 


c d 0 corrected for friction torque, auxiliary 
mechanisms, and so forth 

coefficients in power series for c<i 0 as a 
function of angle of attack 

( c d 0 = s o + + s 2 a r 2 ) 

or 

( c> 1 = 8^ 1 + a + 8 a 
l d Q 0 lr 2 r J 

increment in c^ o to account for friction 
torque, and so forth 

(^do = Cd Q ' " c d Q = A8 0 = 8 0 ' - & 0 ) 
slope of lift curve, per radian 


a 



k 
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T 

Q 


Eotor Aerodynamic Characteristics 
rotor thrust, pounds 
rotor torque, pound- feet 

Q = — A 

5 t(a3 to E)2/ 

Botor-Blade Motion 

flapping angle, radians unless otherwise noted 


k 


a l ,a ’2 , '^l , ^2 




B 1> B 2> B 3> B V B 5> B 6 ,B 7 

m^,a,a> 

Subscripts: 

0 

f 


Miscellaneous 

decay constant in assumed variation of 
induced velocity with time 

constants in expression for 9 as a function 
of time; A]_ and Aq in radians unless 

otherwise noted 

coefficients in differential equation 
for P(t) 

coefficients in complete solution for P(t), 
radians unless otherwise noted 

decay constants in expression for 0(t) 


initial value, for hovering 

final value, for steady autorotative vertical 
descent 


a 


average value 


(c > a - 1[( ) 0 ♦ ( ) J) 


At 


value after time interval of At seconds 
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ANALYSIS 


It is assumed that the initial steady state of hovering flight is 
entirely determined, and the physical quantities of this hovering 
condition are denoted hy the subscript 0. In the final state of 
steady autorotative descent, the subscript f is used to denote the 
physical quantities. 

The simultaneous differential equations governing the transition 
period are six in number: 

(1) For the vertical acceleration of the craft, with the rotor 
thrust a function of ft, V, v, P, P, 6 , and so forth. 

/ 

(2) For the angular acceleration of the rotor, with the rotor 
torque a function of ft, V, v, and so forth. 

(3) For the flapping motion of the blades. 

(4) For the hunting motion of the blades. 

(5) For the induced velocity v as a function of T, Y, Y, and 
so forth ( unsteady flow) . 

(6) For the angle of incidence 9 as a function of the flapping 
and lag angles. 

There is as yet neither theory nor empirical data on which to base 
the fifth of the above equations. For convenience, it is assumed that 
the induced velocity v varies with time after power failure according to 

-Irb 

v = Vf + (v Q - v f )e (1) 

where the value of the decay coefficient k is arbitrary and may be 
assigned different values in order to investigate its importance. It 
may be anticipated that the effect of the approximation represented by 
equation (l) is small, and that the. value of k chosen is not critical, 
when the initial and final values of induced velocity are not widely 
different, as in the usual case. The induced velocity is considered to 
be constant over the rotor disk. While exponential variation of 
induced velocity (equation (l)) represents a convenient assumption, it 
is not unlikely that the actual variation would be somewhat different. 
Yerification of the supposition that the effect is small, .must, failing 
a theoretical approach, await experimental evidence . 
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It is assumed, that the angle of incidence varies with time 
exponentially, and according to 


Q — Q ^ + 


At e 


-a^t 


Ar>e 


-a^t 


( 2 ) 


where A^, A^, a^, and a^ are arbitrary, except that 


h + A 2 ~ 9 0 ~ 9 t 


(2a) 


The changes in incidence due to flapping or hunting are thus 
combined with whatever changes the pilot may make manually. By suitably 
choosing the constants, A^, A^, a^, and a^, rapid or slow pitch reduc- 

tions after power failure may be investigated. Thus, the necessity of 
solving the fourth, fifth, and sixth of the basic equations is obviated, 
by the assumptions represented by equations (l) and (2), and the problem 
simplified enormously. 

It may further be anticipated that the variation of 0 with time 
will be of minor importance both to the designer and in its effect on 
the variations of V and ft with time. For the analysis of flapping 
angle it is therefore assumed that, in the transient period, ft is a 
constant and has the value of the average between ft Q and ft^ . This 

average value is denoted by ft a . 

In the following analysis only untwisted and untapered blades are 
considered. The equations developed on this basis can probably be 
applied to any blades, with fair accuracy, by using equivalent chord 
and angle of incidence. 


Solution for .(3 

The equation for the vertical acceleration of the helicopter may 
then be written as 


P bacB 3 fl a g / 7 , T _ J_ + e 
2 W \2ftgR " 3ft a 3 
g 


V = g - 


(3) 
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and the equation for the transient flapping motion is 



*)- 


^ h 

p acB 0 


a /V 




3fl a R 


P 

w: 


) 




Solving equation (4) for Y, differentiating for Y , and 
substituting in equation ( 3 ) gives 


P + "bq'p + b 2 $ + bjP = 


pacE^o / 

+ — (- T 



pahcB^ftg 2 

+ 4a//g 7 


(5) 


where 


b pacR W b , ^1 N 

1 2 I 1 \I 2 BW / gy 


pacE^fl a / 2 pabcE^jA 
*2 = 6 ^ \ a + 40 W/g / 


h 3 = 


pabcE 2 D a ^ 


W/g 


gpacE^o 


^ = 6 l 


The transient solution of equation ( 5 ) is of the form 


Eh t at 

-I s + 9 h 


P = B, e ^ + e cos cot + B^ sin out^ 


(6) 
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where the roots of 


m 3 + 


b-jm 


+ bgin + b^ = 0 


(7) 


are 


m = 

m = a ± 



(7a) 


where 


j-'Ti 


The steady-state solutions have the same form as the forcing 
functions, and together must satisfy equation ( 5 ) • Thus the steady - 
state solutions are : 


. -kt -ant -aot 

P - B^ + B,_e + Bge 1 + B^e ^ 


(8) 


where, by substitution into equation ( 5 ) and equating coefficients of 
identical terms, 


B 4 = P f 


(8a) 


b 5 = 


jacR^ n_ 

~&U 


L ( T f - T o) 


(* 3 -h 


k + “b^k - 


) 


(8b) 


' P acR 3 fA(IV A i a i " 48w/g A 1 


pabcR^fi 


= 


6 1 


1 / Vl 3 ” ^l 9 -! 2 + *2 a l ‘ ‘ b 3y 


(8c) 
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B 7 = 


pacB 3 fl\ - 48w/g 


pabcE 3 ft 2 
a 


2 " \ a 2^ + h 2 a 2 ' ^ 3 , 


The initial conditions, which must he satisfied by the complete 
solution for p(t), are that, at t = 0, 


P = P 0 


P = 0n = 0 


P = Pn = 


pacE^o^ A"q V 0 . 9 0 


2 H I - p n fl‘ 

3fi 0 E 4 J 0 0 


Therefore, from equations ( 9 ) , 


B 1 + B 2 + B 4 + B 5 + B 6 + B 7 = P 0 


nijB 1 + oB 2 + - kB<_ - a^Bg ~ a 2 B 7 = ® 


(10a) 


2 B 1 + (a 2 - cd 2 )b 2 + 2axfB^ + k 2 ^ + a 1 2 Bg + e^By = P Q tlO 13 ) 


”1^1 


Equations (10), (10a), and (10b) can he solved for B 1 , B 2 , 
and B^, when B^ through By have been found by equations (8). The 
complete solution, then, is 


/x m-, t at / 

P(t) = B^e ■ + e / 


B„ cos aot + B_, sin cut ) + 


-kt -a-it -a^t 

B. + B,-e + Bge + Bye 
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Solution for ft and Y 


The remaining two variables are V(t) and ft(t) . 
be solved are the thrust and torque equations: 


The equations to 


g - 


>abcR3 

2W/g 


<r 


- v A 0 ‘ 

— — ft - “ft + -r0 

2R 3 3 


) 


( 12 ) 


and 


ft = 


pcR^ft 2 

X 2 f x 

. if? 

r f \ f 2 f 3 + 5 0 ' 

. 2ppft 

ro 

H 

H 

2 

3 \ 

3 * f l f 2j - 4 

h 


(13) 


where 


f-^ = a - &2 

f 2 = 0 - P/ft 

f 3 = a p/ft + + f 2 &2 


and p(t) is known approximately from equation (11). 

Equations (12) and (13) can be solved by a tabular, step-by-step 
process as outlined below. 

A sample calculation has been performed (fig. 1), by the step-by- 
step process, comparing the results obtained by including and neglecting 
the p terms in equations (12) and (13) • Reference to the figure will 
show that the effect of flapping is negligible on the variations of ft(t) 
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and Y(t). The flapping terms are therefore neglected in the following 
development of the step-by- step method of solving equations ( 12 ) 
and ( 13 ) . In effect, then, only rigid blades are considered, although 
the effects of blade flapping would be negligible. 

Equations (12) and ( 13 ) then become 


t = g - 


pabcR^/Y - 7 

2W/g \ 2R 


-I- 




( 12 a) 


and 


a = 


pcR^ 

8, ” 

21 

1 1 

3R 2 


(V - - -V-t - ^ 


ft(Y - v) - 


5 0 ' + & 1 0 + 6 2 0 2 
_ 


2 

a 


(13 a) 


By differentiating equation (13a), 


a = 

* 21, * 


(a - O . &, - ©fa - b 0 ) r ... 

2 ~^(Y - T)(Y - v) - -i y SZ[fl( Y - v) + fl(Y - v)j + 


j(a - 25 2 ) 


0( a - 25 J S ' + 8 9 + 8 0 • 

_ — l a ( y - v) - -2 + 205^ 2 
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At zero time (the instant of power failure) all the quantities in 
the above equations (12a), (13a), and (1^-) are known. The _step-by- step 
process is started by computing from these the values of Yq, Qq, 

and 2^ . After a short interval of time, say At seconds, the new value 

of rotor angular velocity o may he found approximately from the 

At 

ahhreviated Taylor series : 


fi At “• n 0 


(At)fl 0 


mi: 


ft i 


(15) 


It is assumed that the change in T in the period At is the 
time At times the average rate of change of T in the time 
interval At. Thus, 


V At - V C 


At ft 

2 VO 


+ T 


At/ 


(16) 


From equation (12a), 



-3 

pahcR J 

^At i 

2W/g 

2B \ 


’^At " v At) + 3 ^At 


(17) 


Solving equations (16) and (17) simultaneously for T 
being tvt and T At^ 


(the unknowns 



g - 


pabcR- 

2W/g 


n At 

2 R 


(vq + ftp - T A t) + 


9 


Ato ^ 
3 


1 + At 

8w/g At 


(18) 


whence 



may he found from equation (16) . 
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Nov 
value of" 


%t 

%t 


may "be computed from equation (13a) and the approximate 
from equation (15) checked hy 


^At = + Tr(% + ^At) 


(19) 


If this does not check quite closely the value first obtained from 
equation (15)# and used in equation (l8), a trial -and -error process can 
be used to find In this case, a value of ^ vould be assumed, 

in place of equation (15)# used in equation (18), and then checked 
by equation (19)* Good agreement should be achieved before the process 
is repeated for the next interval. Experience has shown that 


for At « 1 
2 


second, the first estimate furnished by equation (15) is 


usually checked quite closely (within 0 .05 radian/sec) by equation (19) • 


This process is repeated for as many intervals of time as are 
necessary for ft and Y to reach their final steady values in steady 
autorotation . 


A running check can be kept on the approach to the instability 
caused by excessive blade stalling (see reference l) by computing the 
value of the inflow ratio X at every interval of time, and spotting X 
and 0 on curves of Q against X and 0, in which the effect of 
blade stall has been included, as in figure 3 of reference 1. The value 
of X is 'given by 


X = 


( I Z l) 
ftR 


( 20 ) 


Long before the "second, unstable trim point" (described in 
reference 1) is approached, however, it would be necessary to use 
modified equations (12a) and (13a), altered to include the effects of 
blade stalling. In this case, the trick represented by equation (15) 
to improve the convergence of the step-by-step process would be exces- 
sively complex to apply. The trial- and- error process necessary to 
carry out the numerical integration of the new equations (12a) and (13a) 
would be arduous at each interval of time, and, for adequate accuracy, 
he time intervals would have to be much more frequent . 
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A SMFLE DESIGN STUDIED 

Using the methods described in the preceding section, the tran- 
sitions from hovering to steady autorotation have been computed for a 
helicopter with the following physical characteristics: 

¥ = 2700 pounds c = 1.25 feet 

b = 3 a = 5 06 per radian^ 

B = 20 feet c^' = O.OO87 - 0 00216a,, + 0 .400a r 2 

1^ = 100, 200, 4-00 slug-feet 2 

This sample design is the same as was considered in reference 1, for 
which the steady autorotative characteristics were computed. 

The initial hovering state considered is 

Yq = 0 Vq - 21.25 feet per second 

= 25 ol radians per second - 7 -25° = O v .1265 radian 

For transitions in which the pitch is reduced after power failure, 
the final state considered is 

y = 31.2 feet per second v_ = 25*1 feet per second 

I X 

= 20.8 radians per second 0^ = 4° = O.0698 radian 

For transitions in which the pitch is not reduced, the final steady 
state is taken as 

Y f = 30«5 feet per second v f = 26.4 feet per second 

= 16.5 radians per second = 7 °25° =0.1265 radian 
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The various ways in which the pitch is changed, for this investi- 
gation, are illustrated in figure 2. The different curves correspond 
to different values of the constants A^, A g , a^, and (equation ( 2 )). 

The variations of induced velocity with time, given "by equation (l), 
are illustrated in figure 3, for k = 1, and, in the cases of pitch 
reduction after power failure, for k = 


RESULTS AND DISCUSSION 
Prohahle Effect of Assuming the Variation of 
Induced Velocity with Time 


A typical variation of descending velocity V and rotor angular 
velocity 0, with time after power failure is illustrated in figure 4 . 
These calculations were made hy the step-by-step process given in the 
section entitled "Solution for Q and V," neglecting flapping. Two 
calculations were made, for k = 1.0 and oo^ in order to evaluate the 
importance of k. 

Since, in the usual case, the initial and final values for the 
induced velocity are not very different, it seems intuitively evident 
that the exact variation of v(t) ought not to affect the variations 
of V(t) and ft(t) appreciably . This is confirmed hy figure k } 
where V(t) and 0( t) are compared for the widely different values 
of k, representing the variations .of v(t) shown in figure 3 • 

It may he observed that, if the induced velocity varied in such a 
manner that its value was not always between the initial and final 
values, then it might be important to consider the actual variation „ 

In that case, v(t) would have to be determined from experiment. It 
should be remarked that, in the solution for t) and V(t) just 
referred to, any suitable variation of v(t) may be assumed. The 
exponential variation was convenient in the analytical solution 
for 0(t) . 


Interpretation of the Different Rates of Pitch 
Change Considered 

The extremes of the types of pitch change considered are the upper- 
most and lowermost curves of figure 2, representing no pitch change and 
instantaneous pitch reduction at the instant of power failure . The case 
of no pitch change could occur with rigid blades (or hinged blades so 



l6 


NACA TN 1907 


articulated that a change in lag angle would not change the incidence) 
if the pilot failed to alter the pitch manually. The instantaneous 
pitch change could "be approached by a rotor having articulated blades 
with no damping about the drag hinge. In this case, the change in lag 
angle would be very rapid, and any attendant change of incidence could 
be considered practically instantaneous . The intermediate, exponential 
types of pitch change could each occur, for articulated blades, with an 
infinite number of combinations of pilot reaction and degrees of damping 
about the drag hinge. The more gradual changes are associated with 
slower pilot reaction and greater damping about the drag hinge . For 
rigid blades, the curves represent only different pilot reactions. 


Variations of Flapping Angle P with Time After 
Power Failure 

The influence of the different rates of pitch change, and of blade 
moment of inertia, on the variation of flapping angle p with time is 

given in figure 5 . These calculations were made by the method of 
equation (11), using the induced-velocity decay coefficient k equal 
to unity. 

The variations of {B ( t ) shown in figure 5 are of passing interest 
only, since the effects of blade flapping are of no importance in the 
variations of V(t) and ft(t) , In spite of the linearizing assumption 
in the solution for P(t) that ft is constant, figure 5 probably 

gives a. good indication of stop- settings required on the flapping hinge 
to allow the blades complete freedom in this maneuver. The final values 
of p are different for the case of no pitch change and the cases of 
pitch reduction because the final values of ft are not the same . 


Effect of Blade Flapping on V(t) and ft(t) 

Comparison (fig. l) of two computations of Y(t) . and ft(t), for 
instantaneous pitch change, shows that the effect of blade flapping is 
negligible. The case of instantaneous pitch change was chosen for this 
comparison, since the variations of p are greatest for this 

case (fig. 5) . If the flapping hinge is so directed that a change 
in p causes a change in incidence, then in equation (2) for 9( t) 
the influence of flapping must be considered. The value of k was 
taken to be unity . 
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Influence of Eate of Pitch Eeduction and Blade 


Moment of Inertia on Y(t) and 2( t) 


The influence of rate of pitch change and Blade moment of inertia 
on the variation of flapping angle with time is shown in figure 5; the 
effect of blade moment of inertia on Y(t) and Q( t) is shown in 
figure 6. It may be considered that the average time required to 
effect the transition for the cases considered is about 6 seconds. 

With regard to avoiding the instability caused by blade stalling at 
high values of X (see the study of stability of autorotation, 
reference 1), it is desirable to keep X to a minimum throughout the 
maneuver . 


In order to evaluate the effects of rate of pitch reduction and 
moment of inertia in this regard, figures 7 and 8 are presented. The 
solid lines represent boundaries determined from reference 1, figure 
For any value of 

value of X 


0, the point on the boundary was determined by the 

at which the difference in — — & due to negLecting blade 

a 


stalling was 0 .0015 • Thus the boundaries shown indicate, for any value 
of 0, an arbitrary limit for X, above which blade stalling should be 
accounted for in equations (12a) and (13a). The broken curves are the 
loci of combinations of X and 0 computed during the maneuvers 
considered in the study of the sample design. As can be seen from the 
figures, all the cases computed lie within the range of X and 0 
where blade stalling can be neglected, as has been done in the ANALYSIS . 


From figure 7> the advantage of a rapid pitch reduction is apparent, 
for keeping X to a minimum throughout the maneuver, and hence in 
avoiding the instability of excessive blade stalling. From figure 8, a 
similar advantage for large blade moment of inertia is to be noted. 

For larger initial blade angles in hovering than considered here, 
the maneuver can only be investigated by accounting for blade stalling, 
and then only with grave reservations because of the assumption of 
constant induced velocity, which assumption was shown in reference 1 
to be severe in cases where blade stalling is important. It can, however, 
be anticipated that, in this case, the importance of rapid rate of pitch 
reduction and large blade moment of inertia would be greatly magnified. 


Descending Telocity against Altitude Lost 

Figures 9 and 10 show descending velocity against altitude lost 
during the maneuvers investigated. Altitude lost is obtained by 
graphical integration of the curves of Y(t), figures 5 and 6. The 
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average drop required to effect the transition, for the cases considered., 
is seen to he about 120 feet . 

The effect of rate of pitch change on descending velocity at a 
given amount of altitude lost, from figure 9, is apparently small. This 
conclusion is obviously limited to cases -where the hovering incidence is 
within the range for steady autorotation. If. 0 Q were greater than the 

maximum for steady autorotation (reference 1, fig. 3), then obviously 
the pitch would have to be reduced rapidly to avoid the instability 
caused by excessive blade stalling. 

Prom figure 10, the advantage of a large blade moment of inertia, in 
preserving a minimum descending velocity for a given altitude lost, iB 
readily apparent. 


CONCLUSIONS 


Analytical methods are given for treating the transient motion of 
a helicopter between the hovering state and the steady autorotative 
state of vertical descent. Prom the standpoint of avoiding blade 
stalling in the transition, it iB desirable to reduce pitch rapidly and 
to have a large blade moment of inertia. 

The effects of hinging the blades, for a given net rate of pitch 
reduction after power failure, appear to be negligible. 

The average time to effect the transition from hovering to steady 
autorotation is about 6 seconds, for the cases investigated. The 
corresponding average altitude lost in the transition is about 120 feet. 


Princeton University 

Princeton, N. J., May 12 , 19^8 
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Type of pitch change 

• Instantaneous (flapping considered) 

— Instantaneous 

— Exponential (moderate) 

— Exponential (slow) 

—No pitch change 


Time after power failure, sec 


Time after power failure, sec 


Figure 1.- Effect of rate of pitch reduction on the variations of descending 
velocity and rotor angular velocity with time after power failure. Flapping 

neglected except as noted. I 1 = 200 slug-feet^. 
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Time after power failure, sec 
Figure 2.- Types of pitch change considered. 
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0 2 4 6 8 10 

Time after power failure, sec 


•Figure 3.- Assumed variations of induced velocity with time after power 

failure. 






Rotor angular velocity, a , 

radians/sec Descending velocity, V, ft/sec 


NACA TN 1907 



Time after power failure, sec 



Time after power failure, sec 


Figure 4.- Effect of two different assumptions for induced velocity against 
time on the variations of descending velocity and rotor angular velocity 

with time after power failure, v = v f + (v 0 - v f )e (See fig. 3.) 
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2*— J 1 

0 2 

Time after power failure, sec 


Figure 5.- Effects of rate of pitch reduction and blade moment of inertia 
on the variation of flapping angle with time after power failure. 
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Instantaneous 

• Exponential (moderate) 

Exponential (slow) 

No pitch change 


Blade stalling should be 
considered above 
/, this boundary 
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Figure 8.- Effect of blade moment of inertia on T\ against 0 throughout the 
transition maneuver. Slow exponential pitch change. 




Descending velocity, V, ft/sec 
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Altitude lost, ft 


Figure 9.- Effect of rate of pitch reduction on descending velocity against 
altitude lost. 1^ = 200 slug-feet^. 



Altitude lost, ft 


Figure 10.- Effect of blade moment of inertia on descending velocity against 
altitude lost. Slow exponential pitch change. 






